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^ ■ Abstract 



We study the conservation laws of both the classical and the quantumme- 



! chanical continuum 1/r^ type systems. For the classical case, we introduce 

Q i new integrals of motion along the recent ideas of Shastry and Sutherland 

, (SS), supplementing the usual integrals of motion constructed much earlier 

, by Moser. We show by explicit construction that one set of integrals can be 

related algebraically to the other. The difference of these two sets of integrals 
then gives rise to yet another complete set of integrals of motion. For the 
quantum case, we first need to resum the integrals proposed by Calogero, 
Marchioro and Ragnisco. We give a diagrammatic construction scheme for 
these new integrals, which are the quantum analogues of the classical traces. 
Again we show that there is a relationship between these new integrals and 
the quantum integrals of SS by explicit construction. Finally, we go to the 
asymptotic or low-density limit and derive recursion relations of the two sets 
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of asymptotic integrals. 
03.65.Db 



Typeset using REVT^ 



I. INTRODUCTION 



The integrability of both the classical and the quantum one-dimensional problem of 
particles interacting via the two-body potentials Vo(x) = g^jx^, V^(a;) = 5'^$^ sin~^[$x] and 
V/j(a;) = (7^$^ sinh~^[$x] has been shown more than two decades ago by Moser (for the 
classical problem) and Calogero, Marchioro and Ragnisco ||CMR|| (for the quantum problem), 
both groups exploiting a technique due to Lax [||. These early results have been reviewed, 
extended and collected nicely both for the classical and the quantum cases by Olshanetsky 
and Perelomov in Ref. 

For the classical systems, integrability restricts the motion in terms of action-angle vari- 
ables onto a torus in phase space. However, for the quantum case, integrability leads to 
solvability only for those special cases which support scattering, i.e., systems which fly apart 
when the walls of the box are removed. In these cases, integrability implies conservation of 
individual momenta and thus the wave function is given asymptotically by Bethe's Ansatz. 
For the above interaction potentials, Sutherland |^ has exploited this fact to determine the 
properties of the quantum systems in the thermodynamic limit. 

Recently, Shastry and Sutherland |^S|| have given an independent proof of integrability 
of the quantum many-body problem and constructed new integrals of motion. However, for 
any finite number of particles iV, we know that in principle we have exactly conserved 
quantities. Therefore we expect the new integrals of motion to be related to the integrals 
constructed by CMR. It is the aim of the present work to elucidate some of the features of the 
new integrals of motion and to show their relation to the integrals of CMR. We emphasize 
that this new proof of integrability has also made possible the application of the ideas of 
the asymptotic Bethe Ansatz to the models with quantum exchange 0. 

In section || we show that the new construction of SS gives integrals of motion also for 
the classical problem. We next explicitly calculate these new integrals up to n = 4 and 
compare them to the integrals of CMR. This then gives rise to yet another set of integrals 
Kn- Section |I| is devoted to a comparison of the two series of integrals of motion for the 
quantum case. The integrals given by CMR are not extensive quantities and we need to 
resum them via an application of the linked cluster theorem. In section |V], we take the 
asymptotic or low-density limit of the problem and section summarizes and discusses our 
results. 



II. THE CLASSICAL CASE 

The Hamiltonian of primary interest for our present work is given as 

a = Y.p1 + - l)$'Esinh-'[$(x, - X,)]. (1) 

i ij 

The interaction term reduces to Vq in the limit of high-densities (or $ — 0) and the trigono- 
metric interaction Vt is just the analytic continuation of $ ^ i$. Here and in the following, 
we will use the primed sum Yl! to indicate that the summation runs over unequal indices 
only. 
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A. Moser's invariants 



Let us briefly recall the method of |jl|,@]: We introduce the Lax pair L, M, 



Ljk = PjSjk + i{l - 5jfc)v^(^ - l)*^ coth[$(xj - Xk)], (2) 

6jk J2 smh-2[$(xj - xi)] + (1 - 6jk) smh-^[^{xj - Xk)] . (3) 



M,fc = 2^A(A-l)$^ 

L I J 

The classical equations of motion then imply the matrix equation 

^ = {L,H} = t[ML-LM], (4) 

where we define the Poisson brackets as {F,G} = Ylf=i ~ §^§§'- The time evo- 

lution of L consequently is an 

isospectral deformation, L{t) = exp[i /q M(r)(ir]L(0) exp[— Jq M{T)dT]^ and the integrals 
of motion are simply given as the traces 

r„ = TrL"(t) = TrL". (5) 

We also need to show that the T^'s are in involution, e.g., {T„,Tm} = 0. Using the 
Jacobi relation for Poisson brackets, we see that 

{H, {T„, T„}} = {T„, {H, T„}} - {T„, {H, T„}}, (6) 

and thus {T„, T^} is also an integral of motion. But, allowing the system the evolve in time, 
all particles scatter, the Lax matrix itself evolves into 

T ■ — h- 



L'^L°^ = \ Ljk = +t^XiX - 1), j > k, (7) 
Ljk = -i\J\i\ - 1), j < k. 

and so the coordinate dependence vanishes. Thus the Poisson bracket {T„,Tm} evaluates to 
zero. We remark that it is this procedure that we use to prove involution for all the integrals 
constructed in the following chapters. 



Let us define ajk = yA(A — 1)$ coth[$(a;j — Xk)] and ajj = 0. Then a direct calculation 
of the integrals of motion up to n = 4 gives 

Ti = ^p, = P, (8a) 

i 

T2 = T.Pl + t<^l (8b) 

i ij 

i ij 
I 

T, = J2pt + ^J2 4(P? + m + P') + (8d) 
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Using 



$^A(A - 1) [smh-2[$(xj - Xfc)] + ij , we see that = H+<I>^X{X - 
Note that due to the antisymmetry ctjj = — Oji, only even powers of a — and thus integer 
powers of A — will appear in all these expressions. 
Let us now define the classical down-boost |H1 



X 



N 



We then find easily that 



{X,T„} = nTn-i- 

Further, Jacobi's identity gives 

{X, {T„, Tm}} = {n- l){T„_i, T^} + (m - l){T^_i, T^. 



(9) 



(10) 



(11) 



As a particular case, suppose n = 2, son — 1 = 1 and T^^i = P. Then by translation 
invariance {P, T„} = 0, so we conclude that if {if, T„} = 0, then {if, T„_i} = 0. In 
particular, {ii, T/v} = implies that all T„ are integrals. Finally, we may construct all 
integrals of motion from T/v by repeatedly using the boost X in the representation 



X 



(12) 



B. Shastry's invariants 

In Ref. I^Sf , Shastry and Sutherland provide a set of integrals of motion for the quantum 
problem. However, mimicking their arguments, we can straightforwardly construct integrals 
of motion for the classical case, too. Let us introduce the singular matrix Aj^ = 1 for all 
z, j and the vector rjj = 1 for all j. Then we define integrals of motion s.t. 

J„ = Tr[L"(t)A] = r^tL-(t)ry = ^ L,,,,L,,,, ■ ■ ■ L,„_,„ L,„,„^„ (13) 

il,i2,...,in + l 

with the Lax matrix L given as before. We then have 

^ = 1 {Tr [exp[zMt]i:"(t) exp[-iMt]A]} (14) 

= 2Tr[Mi:"(t)A-L"(t)MA] (15) 
= I [Tr[L"(t)AM] - Tr[L"(t)MA]] (16) 
= 0, (17) 

since MA = AM = as shown in SS. Involution for these integrals of motion is proven by 
the same asymptotic argument as before. A direct calculation of the conserved quantities of 
SS up to n = 4 gives: 
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i 

I I 



ijk 



(18a) 
(18b) 
(18c) 



ijk 



Ji = Ypt + '^Y (^Ij Pi + PiPj +P] + Tra"^ + ^ <yij<yjk<ykiOiim 

i ij ij^j^kj^l^m^i 

I 

- aijajk [p- + + Pfc + PiPj + PjPk + PkPi 

ijk 



(18d) 



Again the Hamiltonian can be found in the n = 2 term, J2 = H + $^A(A — 1)A^(A^^ — l)/3 
and again only even powers of a appear in the expressions of the J„'s. 

The action of the down-boost on these new integrals of motion is as in Eq. (pUj), e.g., 
{X, Jn} = nJn-i- Much more useful is the up-boost Y which we define as 



(19) 



in analogy with the up-boost operator J2n ^S^Sn+i in the Heisenberg model. Unfortunately, 
this up-boost only works, if we restrict ourselves to the potential Vq such that afj = A(A — 
1) / (xj — Xj)"^. In this case, we find by explicit construction that {Y, Jn} = {n + 1) J^i+i- The 
Jacobi relation { J„, {Y, J„}} = {Y, { J„, J„}} - { J„, { J„, Y}} now gives 



{n + l){Jm, Jn+l} = {Jm, Jn}} - (m + 1){ J„, Jm+l} ■ 



(20) 



Thus, if {Jm, Jn} = and {Jm+i, Jn} = 0, we also have {Jm, Jn+i} = 0. We emphasize that 
the up-boost (^) seems to work only for the special potential Vq. 



C. Relation between invariants 

We can again use the Jacobi relation to show that the Poisson bracket {T„, J^} is an 
integral of motion which in the asymptotic limits evaluates to zero. The difference between 
the integrals of motion of Moser and SS then gives rise to yet another set of constants, 

Kn = Jn ~ Tn = Li-^i^Li^i^ ■ ■ ■ Li^_j^i^Li^i^_^-^^. (21) 

Various terms in the J„'s can be simplified with the help of the coth-rule, 

aijajk + aijOki + ajuau = -$^A(A - 1), (22) 

and hence we find 
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K3 = $2A(A 
= $2a(A 



1)N{N -1){N -2)/3, 
l){N-l)iN-2)P, 

l){N-2) [(A^ - 2)T2 + P^l + [<I>2A(A-1)(A^-1)(A^- 2)1^9 



(23a) 
(23b) 
(23c) 

(23d) 



Note that is the first term that is not a simple constant, and in order to make the K^s a 
complete set of integrals of motion, we may simply use K^+i and -R'Ar+2- Thus we conclude 
that by construction, we can express Shastry's integrals of motion in terms of Moser's and 
vice versa. We emphasize that this relationship is not linear, but only algebraic as seen from 
the existence of the term in K4. 

Taking the limit $ — > 0, we see that the K^s are zero. Thus only for the simplest case 
of the Calogero potential Vq{x) = /x^ do we find that the Moser set of integrals of motion 
is identical to the set of SS. 



III. THE QUANTUM CASE 

In the quantum case, the elements of the Lax L and M matrices become operators 
themselves, i.e., the momentum operator is pj = —id/dxj and we have the commutation 
relation [xj,pk] = iSjk- Since operator elements do not necessarily commute, we always 
mean an ordered product of elements when we multiply matrices in the following. 



A. Calogero's invariants 

The early work of Calogero, Marchioro and Ragnisco ||CMR|| quantised the classical Lax 
equation, by antisymmetrizing the right-hand side of Eq. (^). The proof of invariance of 
the traces then does no longer hold. However, CMR also showed that after replacing the 
classical variables with the corresponding quantum mechanical operators, we can define new 
integrals of motion /„ s.t. 



N 



A{P) = det[l -PL] = 1 + Y.i-PT^n 



(24) 



n=l 



CMR then go on to argue that these /„ are conserved, [/„, H] = 0, and in involution, 
[In, Im] = 0. The later result is again proved [0 by use of the asymptotic limit as in the last 
section. A direct calculation of the conserved quantities of CMR up to n = 5 for H yields: 



h = i/f - i 

2^2 



a. 



1 ' 1 ' 

^3 = g E PiPiPk - 2 E 
ijk ijk 



(25a) 
(25b) 
(25c) 
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1 ' 1 ' 

h = -^ J2PiPiPkPi - 4 E ^^PkPi 

ijkl ijkl 
1 ' 1 ' 

~ 4 51 aij^jkakiaii + g H (25d) 

1 ' 1 ' 

^5 = ^ II PiPjPkPlPm - Y2 ^ C^ljPkPlPm 
ijklm ijklm 

1 ' 1 ' 

~ 4 XI OiijO^jkakiaiiPm + g XI ^'ij'^liPm- (25e) 

ijklm ijklm 

Note that the Hamihonian can be found in the term in parenthesis in J2. 

Let us define a quantum down-boost operator analogous to the classical boost 0] . With 
X = Y^jLi Xj as before, we then find 

[X,I^]=t{N-m + l)Irn-i. (26) 

Using Jacobi's identity for commutators, we can easily show that as previously, [H, In] = 
implies [H, = and thus [H, /jv] = implies all J„ are integrals. A particularly nice 
result is to write In = det L, treat the momenta pj as classical c-numbers since there are no 
ordering ambiguities, and use the representation 

TV g 

to generate all /„ in the quantum case. 

Of special importance in the following will be that as in the classical invariants by Moser, 
a will only appear in even powers in the J„'s. Therefore, A will occur with integer powers 
only and terms such as [A( A — 1)]^/^ do not exist. 



B. Shastry's invariants 



In Ref. |pS|| , Shastry and Sutherland provide a proof of integrability in the quantum case 
via an entirely different method: The Hamiltonian H is given as before but the Lax matrices 
now read 



Ljk = PjSjk + ^(1 - Sjk)\^ coth[$(xj - Xk)] 



= Pj^jk + ^(1 - ^jk)Xjh 



5jfc^sinh 



Xl 



6jk) sinh - Xk) 



(28) 
(29) 

(30) 



with Xii = 0. SS define their conserved quantum integrals of motion as in Eq. (p!8D , e.g., 
Jn = vi\L^^)^ti. The new Lax matrices obey the ordered Lax equation 



and we may easily prove invariance via 



(31) 
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[J„, H] = v^liL^^T, H]r] = r/t [m^^(L^^)" - (L^^)"M^^1 7^ = 0, 



(32) 



since as before Vj'^M^^ = M^^rj = 0. A direct calculation of the conserved quantities of SS 
up to n = 4 yields: 



i 

I I 



ijk 



•^3 = ^ + 3 J2ixl + Xij)Pi - XijXjkiPi + Pj + Pk), 



(33a) 
(33b) 
(33c) 



ijk 



Ji = Y,Pi+'^ J2ixl + x'ij) [p1 + PiPj + PjJ + XijXjkXklXlr. 
i ij i^j^k^lj^m^i 

I 

+Trx^ - Y XijXjk [pI +p]+pI+ PiPj + PjPk + PkPi 

ijk 

I I I 

+2^ XijPj + 4i ^ XijX'ijPj + i Y XijX'jkiPj - Pk) 

ijk 



Y X^j -'^Y XijXij + 2 ^ XijX'jk 

ij ijk 



^ X Xti) ~l~ ^ y XijXjk ~^ ^ y XijXjk ~^ 2 ^ ^ XijXij 



ijfc 



- H [xiiXj-fcXfci + "^XijXjkX'ki] + Y XijXjkX'ki 

ijkl ijk 



(33d) 



The derivative x'jk defined by the commutator [pj,X^jk] = ~'^X%^^^ ■ See the appendix for 
an explicit list of derivatives. 

Using x'ij = — sinh~^[$(xj — Xj)], we see that just as in the classical case, J2 contains 
the Hamiltonian, i.e., J2 = H+^'^X'^N{N'^ — l)/3. However, the interaction strength A(A — 1) 
in the Hamiltonian could only be obtained with the modified form of the Lax matrix 
Also, the A dependence of the constant term in the above equation is different from its 
classical counterpart. We remark that the last terms in Eq. (|33b|) and ( P3c| ) can again be 
written as const, and const, x J2iPi by the coth-rule of Eq. (^). 

The down-boost operator acts as before, e.g., [X, J„] = inJn-i- In case of the potential 
Vq, we may also use the up-boost of Eq. ([191) in operator form as 



Y = Y(^d + Pi^^)/2 + Y(^^ + ^.)4/2. 



(34) 



Then [F, J„] = i{n + 1) J„+i and we again have from the Jacobi identity 
i{n + l)[Jm, Jn+i] = [Y, [Jm, Jn]] - i{m + l)[J„, Jm+l], 



(35) 
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so if [Jm, Jn] = and [Jm+i, Jn] = 0, this then imphes [Jm, Jn+i] = 0. We remark that an 
operator similar to our up-boost operator Y, which we constructed in analogy to the boost 
in the Heisenberg model, has been found previously by Wadati, Hikami and Ujino in the 
context of an investigation of the systems with algebraic potential Vq 0. 

Finally, we note another interesting property of these integrals of motion: Let \E'o denote 
the ground state of the model, then it has been shown in Ref. |§] that J2j Lff'^o = for all 
i = 1, . . . , N. Therefore, we see that 

^Un^O = (36) 

for all n. Thus all the J„'s somehow know about the ground state and subtract the appro- 
priate expectation values, e.g., the ground state expectation value of the Hamiltonian is just 
the above constant ^^X'^N{N'^ - l)/3. 



C. Perturbation theory in the Lax matrices 



Looking at Eq. (pSD, we see that each /„, n > 1 in fact contains various powers of /i. 
Furthermore, in the thermodynamic limit, the J^'s are not extensive quantities. Thus the 
situation seems to be similar to the usual problem of connected and disconnected pieces of 
diagrams encountered in perturbation theory. In brief, CMR's /„ seems to contain discon- 
nected pieces and we hope that by a linked cluster expansion, we can write new integrals of 
motion with connected graphs only. 

Let us be specific: With the help of the fermionic coherent path integral |Tp[, we may 
rewrite the determinant 

A(/5) =det[l-/5L], (37) 



/ 
/ 



Yl dcldca exp [- c* [6jk - /3Ljk]ck] , (38) 

a jk 

Y[ dcldca exp [-P S* i^jk/P)ck - c*LjkCk] , (39) 

a jk 



where c*,C(j,a = 1,...,A^ are Grassmann variables. Note first that we may write this 
expression both for a classical L and a quantum L. The fact that the elements of a quantum 
matrix will not necessarily commute with each other is been taken care of by the Grassmann 
nature of the integration: each momentum pi will only encounter indices j ^ i, otherwise 
the integration measure will have expressions like CjCj or c*c* which are zero. 

When we now include a dummy time dependence for the Grassmann variables, i.e. c^*^ = 
ci*"* (i) ) can write 



A{(3) = f l[dc:{T)dca{T)exp[- f dt{Yc;{t){6,k/f3)ck{t) - c*{t)L,kCk{tm (40) 

a "'0 jk 

= Ao(exp[- f dtC£~c]{t)L,kCk{t))])o. (41) 



dt{ 

jk 

where the average is defined as 
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{F{c:iU)clit,)...c,{h)c,{U)...))o 



A, 



^ ''lldc:>{r)dCa'{T)exp[- f dtY,c}{t){l/ f3)c,,{t)] x 



F{cl{U)cl{t^) . . .Cg{h)cn{ti) . . .). (42) 

This is very much hke a path integral description of a many-body partition function Z. We 
further note that the interaction part V = J2jk (^*j(t)LjkCk(t) is just the super Lax operator 
C of SS. 

The perturbation expansion is obtained by expanding Eq. ( ^OD in a power series 



A(/3)/Ao = E / dhdh ...dtn 

(Yl ■■■Y. <(^n)^WnCj„(tn))o, (43) 



- E (44) 

n=0 

and A,„ ~ The last equation is obtained by comparison with Eq. ( ^If ) and 

Ao = 1. Note that /„ = in Eq. (||) for all n > N. E.g. for = 2, we have 
h ~ EfjXj2 1^2? ^nii^i2i2^i3i3(c*iCiiC-2 0^240^3)0 and clearly 23,^3 always take index val- 
ues already covered by {ii, ji,i2, j2}- Thus the bracket () is zero by the Grassmann character 
of the c's. 

Let us now calculate the first few orders of A(/5). With gi being a dummy propagator, 
we find 

Ai = -pY.Lii9i (45a) 

i 

^2 = E (LijLji - LiiLjj) gigj (45b) 

Jl J2 J3 

= ~^/^^E i^ikLjjLki — LikLjiLkj — LijLjkLki 

ijk 

+ LiiLjkLkj + LijLjiLkk ~ LiiLjjLkk) 9i9j9k- (45c) 

Introducing the diagrammatic notation i» — > j = Lij, we can represent these expressions 
by their graphs as in Fig. Note that only in Eq. ( [45(j ) do we need to worry about the 
ordering of the matrix products. If we ignore that ordering for the moment — the classical 
case — we have 

A3 = "■^/^^ E [^^ii^jkLkj — '^LijLjkLki — LiiLjjLkk] gig^gk- (46) 

ijk 

We see that the second term in Eq. (|46|), representing the fully connected diagram, is actually 
just const. X Tr(L^). 
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Let us see what these expressions tell us: (i) We see that in fact the A„'s are just the 
/„'s of CMR with Qi = 1. (ii) We observe that the fully connected diagrams give the traces 
of powers of L in the classical case as was expected from the well-known matrix formula 
Indet A = TrlnA. Thus these connected diagrams are the quantum analogue of the classical 
integrals of motion, (iii) The ordering of the matrix products becomes important for n > 2, 
thus necessitating order labeling of diagrams. 



D. Constructing connected diagrams 



We now want to rewrite the perturbation expansion (^) such that we only use fully 
connected diagrams. And we want to do this such that we can minimise the ordering 
problems coming from the quantum character of the Lax matrix. The basic program is due 
to Thiele and know as the linked cluster theorem. It can be summarised as follows: We 



resum the series (|^) as 



Ai/3) = exp[- Y: 



/3" 



n=l 



n-1] 



-T 



and use it to define the T„'s. Comparing Eq. (^if) and (13) 5 find up to n = 4 
Ti = Ji, 



To 



It - 2/2, 

i ij 
If - 3/1/2 + 3/3, 

i ij 

/f- 4/2/1' + 2/2 + 4/3/1 -4/4, 

Y.Pt + ^J2 (^Ij \v] + ViVj + v1\ + Tra 



(47) 

(48a 
(48b 

(48c 
(48d 

(48e 
(48f 

(48g 



2 ^ Oi^jQij 



2EI 



a. 



AiJ^o^ijO^ijPi- 



(48h) 



Since CMR have already proven [/„, Im] = 0, there is no ordering problem for the /„'s in the 
construction of the T„'s and we furthermore have [T„, T^] = 0. Since T2 is, up to a constant, 
the Hamiltonian, this implies both involution and invariance. 

As expected, we find that each T„ corresponds to the fully connected diagrams of the 



series (|4^). We can now directly use the diagrammatic approach to construct the T„'s. 
However, for a given n, there are (n — 1)! different labeled diagrams and thus different matrix 
orderings. Each diagram itself is an ordered operator expression and it is quite tedious to 
get them into a form as in Eq. (^) with all momenta to the right. As an example, we give 
the diagrams for T4 in Fig. ^. 
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Ignoring matrix and quantum ordering, the resultant expressions for the T„'s are equal 
to the classically invariant traces of Eq. (^. Thus we may hope that due to the special 
form of the quantum Lax matrix L, the matrix product order somehow is unimportant and 
the T„'s are just the quantum traces TrL". Explicitly calculating the quantum traces, we 
find that indeed up to n = 3, we have T„ = TrL". However, for n = 4, the quantum trace 
includes the nonzero term —2Y,[ji^aijajka'i^i- Note that this term has a factor [A(A — 1)]^/^ 
TT| . But as shown in section |111 A| , such a term does not arise in the J^'s and consequently 



also not in the T„'s. Therefore, the T^s are not simply the quantum traces. 

Note that these expressions are again very close to the ones obtained for J„. However, 
as before in the classical case, we see that already for n = 2,3, there are the same constants 
in the J expressions which do not appear in the T expressions. 



E. Relation between invariants 

We again would like to see if we can express Shastry's integrals in terms of our T^'s. As 
we have seen in section [11 C| for the classical case, we expect this relation to be algebraic. 



Fortunately, as shown in the last section, fractional powers of A neither appear in the J„'s 
nor in the T„'s so that no a priori reasons forbid an algebraic relationship in the quantum 
case. 

Furthermore, given two sets of integrals of motion {T„} and {Jn}, we know that com- 
mutators of integrals are themselves integrals of motion, and since asymptotically these 
integrals evaluate to zero, the two sets of integrals can be simultaneously diagonalized. A 
relationship between asymptotic integrals of the form J7„ = y4„[{7^}] can always be found, 
since either set of integrals gives an algebraically complete set of symmetric polynomials 
of increasing degree. Suppose we have such a relationship. Then, the operators Jn and 
^n[{^m}] have the same eigenvalues in the same basis, hence must be the same operator, 
and so there must exist a relationship J„ = y4„[{Tm}] between the operators themselves. 

Replacing aij and Xij by their appropriate definitions, using the explicit form for the 
derivatives as given in the appendix and counting powers of A and p's, we then have 

Ji=Ti, (49a) 
J2 = $2A[3 + A(iV - 2)]N{N - l)/3 + T2, (49b) 
J3 = $'A[3 + X{N - 2)]{N - l)Ti + Ta, (49c) 
J4 = [-75 + 120A - 48A2 + (55 - llOA + 52A^)iV + (-5 + 25A - 18X'^)N^ + 2A^A^^] x 
A^A^(A^ - l)$Vl5 + [2 + A(A^ - 2)] X^^T^ + 

10 - 12A + 11(A - 1)A^ + (2A - 1)N^] A^^Ta + T4. (49d) 

Hence we have succeeded in writing Shastry's quantum integrals in terms of the T„'s which 
in turn are derived from Calogero's quantum integrals for up to n = 4. Again, as in the 
classical case, this relationship is not linear, since we observe the term in Eq. ( |49dD . And 
only if we restrict ourselves to the potential Vo(x) = g"^ /x^ by taking the limit $ — 0, do 
we find that both sets of integrals of motion are identical. 

With the A^-particle ground state as before, we may use Eq. (|36|) and hence relate 
the expectation values of various T„'s. E.g., ^JTs^o = $^A[3 + A(A^ - 2)]N{N - l)/3 and 
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-$2A[3 + A(A^-2)](Ar-l)^jTi^o ~ ^o^'^o = 0. We further note that as in the 
classical case, we may define new non trivial constants of motion Kn = Jn—Tn for $ 7^ 0. We 
then have '^iKn'^Q = — \E'qT„\E'o- Unfortunately, we can not give a simple formula directly 



in terms of the Lax matrices for the construction of the KnS analogous to Eq. (21). 



IV. THE ASYMPTOTIC LIMIT 

In the asymptotic t — 00 limit — equivalent to the low-density limit — the elements of 
the Lax matrix are no longer operators, but numbers. Thus explicit calculations are much 
easier and we hope that we can study the connection between asymptotic Calogero and 
Shastry integrals of motion in more detail than in the last section. 



A. Asymptotics of Calogero's integrals 



The asymptotic form for the Lax matrix L gives a corresponding asymptotic form for 
the Calogero integrals /„ — > X„. We define a generating function of the asymptotic Calogero 
integrals. 



N N 

Y[il- zpj - iXz) + Y{{1 - zpj + iXz) 



and then have 

det[l - zL°°] = V[z, ^A(A-l)]. 
We also define the elementary symmetric functions of the N variables pj as 



N 



Orb] 



l<il<-<jr 



and for convenience Oq = 1, and a_r = 0. Then 

00 

V[z,\] = Y.{-\^z'YaN,2r[l-zp]. 

r=0 

In particular. 

In = aN[p] - A(A - l)aN~2[p] + A^(A - l^aN-ilp] 
and the other integrals can be constructed via 



■ N g 



(iV-j + l)J,_i. 



(50) 



(51) 



(52) 



(53) 



(54) 



(55) 



Since the elementary symmetric functions aj obey the same relationship, this also gives us 
the expression for Ij as a linear combination of The general expression 

for Ij in terms of the ar[p] can be obtained using 
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[i-^Pi = i:(-^r^f^"^/ '■)<■.-», (56) 



j=0 \ -I 

and we find 



^.-i:[-A(A-l)rf^+?' ^]a,.,M- (57) 



B. Some generating functions 

Let us define the quantity 

N 

Z[z, A] = - ^(Pi + ^A)] = V[z, A] - iA;2A/'[^, A]. (58) 

The symmetric part ol Z \sV = {Z[z, \] + Z[z, —X\) /2 as in the previous section, while the 
antisymmetric part (^^[2;, A] — Z[z,—X\)/2 is given by 

00 

iXzN[z, A] = iXz Y,{->^^z^yaN-2r-i[^ - zp]. (59) 

r=0 

For z real, these are the real and imaginary parts of Z, respectively. This expression for Z 
is the standard generating function for the elementary symmetric functions, so that 

N 

Z[z,X]=J2i-^y<^j\P + ^>^]- (60) 

j=0 



Clearly, then, we have Ij = Ke{aj\p + iyX{X — 1)]}. Taking the logarithm of Z, 

N 

In Z[z, X] = J2 - z{p + j + ^A)], (61) 

one advantage of using the generating function Z is clear when one anticipates the thermo- 
dynamic limit. We now consider 

dz 1 - 4Pj + ^A) 

This, however, is the generating function for the symmetric power sums 

N 

br\p] = Epj (63) 

since 



(62) 



r=0 7=1 



ZPj 

and we see that 



(64) 



— In Z[z, X] = J2 z''br+i\p + iX] = P[z, A]. (65) 

r=0 
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C. Asymptotics of Shastry's integrals 



Shastry's integrals also approach an asymptotic form Jj — > J}, and similarly we define 
an asymptotic generating function 

_ U[z, A] 

then, we see that jV'f-z, A] = ^[z, AlPf^;, A], so expanding, we have 

A'' oo oo 

Y^{-zfh^{ak\p + i\]} = A^^(-iy-ij,_iRe{afc_,[p + iA]}. (67) 

k=l k=lj=l 

Equating powers of this gives 

k 

Im{afeb + «A]} - A^(-iy-^J,_iRe{afe_,[p + iA]}. (68) 

i=i 

More explicitly, 

X~Hm{ak\p + iA]} = NRe{ak[p + - JiRc{ak^i[p + iA]} + . . . + {-ifjk- (69) 
This allows us to iterate and find Jk in terms of the other Jj, 

Jk = Jfe_iRe{ai[p + ^A]} - ... - {-lfNRe{ak[p + + (-l)'A-^Im{afc+ib + ^A]}. (70) 

Since 1j = Re{aj[p + this recursion relation also relates the asymptotic forms of 
the CMR and SS integrals. However, this relation is between asymptotic integrals with the 
same parameter, and such integrals do not even commute. 

V. CONCLUSIONS 

Our original intent and hope at taking up the present work was to give a simple connec- 
tion between the integrals of motion of CMR and the recently discovered integrals of SS. In 
fact, we were speculating that due to the special structure of the quantum Lax matrix, we 
would simply find T„ ~ J„. 

Quite the opposite has happened: In the quantum case, though we have succeeded in 
rewriting the /„ integrals of CMR into extensive quantities T„, we have however failed to 
give a simple formula for the connection of these T„'s to the J„ integrals of SS for all but the 
simplest potential Vo{x) = g^/x^. In general, we do find a complicated algebraic relationship 
which gives rise to yet another set of non trivial integrals of motion K^- 

In the classical case, we show that the quantum definition of SS may also be used to 
construct classical integrals of motion. Hence here the situation now is just as in the quantum 
case and again we show that we may reexpress the integrals of SS in terms of the classical 
integrals of Moser. Again, the difference of these two sets of integrals vanishes only for the 
potential Vq and otherwise may be used to define new constants and this time, we can give 
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an explicit formula for the KnS directly in terms on the Lax L matrix. Indeed, we may even 
define a one-parameter family of integrals of motion in the classical case, e.g., 

i?„(<5)=Tr[L"(l + 5A')], (71) 

with A'-j = 1 if i j and zero otherwise. Then i?„ interpolates between Moser's integrals 
-Rn(O) = Tn and the integrals of SS -Rn(l) = Jn- 

Most of the previous formulas are given in terms of a's and x's and are thus valid not 
only for the hyperbolic pair potential Vh, but also for the trigonometric Vt and the rational 
Vo after taking the appropriate limits as mentioned in the beginning of section 0. However, 
we have found an up-boost only for the classical and the quantum many-body system with 
algebraic potential Vo(r) = g'^/r'^. Also, the elliptic potential = 1/sn^, is not included 
in this study: Although the classical integrals of Moser and the quantum integrals of CMR 
are valid for this potential, the proof of SS does no longer hold both for the classical and 
the quantum case. This is due to the fact that the row and column sums of the elliptic 
Lax M matrix are no longer zero. The Ansatz = TrL"(7[{xj}] gives an equation for the 
coordinate dependent matrix g as -^g = Mg — gM. Unfortunately, we have not found a non 
trivial such g such that it gives the SS integrals of motion in the elliptic case. 
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APPENDIX A: DERIVATIVES 

We have defined = coth[$(a;i — Xj)] and so Xij = ^^lij and aij = $-y/A(A — l)7jj. 
Then the derivatives are 

^ = x;=$^A[l-7a (Ala) 

x:;. = -2<|.3A7.,(l-7j) (Alb) 
x::; = -2$^A(l-47j + 37j). (Ale) 



The same relations hold for aij, i.e. a- - = $^ JA(A - a-'- = -2$^WA(A - l)7jj(l 



7|,), anda- = -2<I>YA(A - 1) (1-472+37,^.), and, lastly, 7.', = $[1-7^,], l", = -2$7.,(1- 
7^),and7l^ = -2$(l-47|,+37^). 
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FIGURES 



FIG. 1. Diagrammatic representation of Eq. (^). Each line is labeled by a number indicat 
precedence in the corresponding matrix product. 1 corresponds to the right-most matrix. 

FIG. 2. The 6 diagrams associated with T4. 
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Figure 1, R.A. Romer et al., "Conservation laws in the continuum systems". 
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Figure 2, R.A. Romer et al., "Conservation laws in the continuum systems". 
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